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ABSTRACT 

Energy and diffusive mass transport associated with the thermospheric 
circulation are considered in a self-consistent, though mathematically relatively 
simple form to describe in a three-dimensional two- constituent model magnetic 
storm characteristics in composition (N 2 , O, and He), temperature and mass- 
density. It is shown that during disturbed conditions the latitudinal variations 
of composition and gas temperature T g reflect the local nature of the magnetic 
storm heat input assumed to be primarily confined to the auroral zones. There- 
by T and N increase, He decreases and O remains constant through the 
auroral zones at exospheric heights (due to the superposition of temperature and 
diffusion effects) in agreement with OGO-6 mass spectrometer measurements. 

In contrast, the magnetic storm response in the total mass density is character- 
ized by a strong world-wide component and a relatively insignificant increase 
toward the poles with the density peak occurring between two (poles) and eight 
(equator) hours after the maximum energy input, in substantial agreement with 
satellite drag data. While in situ composition and satellite drag mass-density 
measurements can thus be reconciled it must, however, be emphasized that the 
temperature derivation from the satellite drag data cannot be justified during 


disturbed conditions. 
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I. INTRODUCTION 

Jacchia (1959) and Jacchia et al., (1967) have shown that the thermospheric 
density increases during geomagnetic disturbances with the effect being propor- 
tional to the planetary index A p for large disturbances and to the index K for 
small disturbances. Furthermore it has been observed from satellite drag data 
that the effect, occurs worldwide but is systematically stronger at higher latitudes 
(Roemer, 1967). 

In a recent paper by Taeusch et al. (1971) OGO 6 mass spectrometer data 
were presented which revealed at 400 km: 

(a) a gradual increase in N 2 toward high latitudes (up to a factor of 
10 for A p = 130), 

(b) an increase in the exospheric temperature (from 1000 to 1400°K) at 
high latitudes which was inferred from the variation in N 2 , and 

(c) negligible latitude variations in O. 

Volland and Mayr (1971) treated the magnetic storm response of the ther- 
mosphere in a simple analytical model. Assuming that the energy source is 
confined to auroral zones at ±65° we expanded this input into a series 
of spherical harmonics ( P ) 

Q = 2q n P n =Q 0 (P 0 + 3.66 P 2 + 2.22 P 4 - 2.65P fi ••■) 
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which, as expected, converges slowly. However, due to the importance of ad- 
vective energy exchange, associated with the thermospheric circulation, the ef- 
ficiency for the excitation of density variations decreases with the power of 
about l/n 2 and consequently higher order terms produce only small contributions 
for the density structure. This characteristic thus explained the observed world- 
wide component as well as the gradual increase toward higher latitudes in the 
thermospheric mass density during magnetic storms. 

It was shown by Mayr and Volland (1970, 1972a) that during disturbed condi- 
tions the thermospheric circulation effectively redistributes the minor consti- 
tuents such as to reduce the oxygen concentration at high latitudes and within 
the lower thermosphere. Since this depletion is associated with a thermal ex- 
pansion both effects tend to cancel each other in the upper thermosphere hence ex- 
plaining the negligible response in O at times when the N 2 concentration is ob- 
served to increase (Taeusch et al., 1971). In that paper we restricted ourselves 
to a very simple model by postulating that in a first approximation 

(a) the latitude dependence in the form of P 2 (0) dominates, 

(b) a diffusion model can be adopted in which the momentum feedback from 
the minor to the major constituent is negligible, and 

(c) the effect of the diffusive redistribution can be neglected for the energetics. 

Although none of these simplifying assumptions affected our earlier conclu- 
sions substantially, it will be shown in this paper that by relaxing them a number 
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of important features evolve which further characterize the magnetic storm 
response in the thermosphere. 

II. QUALITATIVE DISCUSSION 

Figure 1 shows a simplified block diagram for the thermosphere dynamics. 

It illustrates the links between the composition and gas temperature on one hand 
and the various physical processes and atmospheric parameters on the other hand. 

Suppose the heat input Q is known as a function of time and space defined with 
height and latitude. Some of the energy is conducted down into the lower atmo- 
sphere thus affecting the temperature distribution and with that the composition. 

However, a second and not less important energy drain goes through the 
global circulation which is significantly affected by ion collisions. In this me- 
chanism which is explicitly neglected in one dimensional models energy is ex- 
changed through advection and this in turn affects the temperature and the 
thermospheric composition. 

In parallel the thermospheric circulation induces diffusion which depends 
to various degrees upon the eddy diffusion coefficient, and this diffusion process 
affects the composition directly. Any variations in the composition are 
associated with variations in the partial and total pressure fields. This in 
turn influences the global circulation and with that advection and with that 
again temperature and composition. So it is apparent that all the processes 
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surrounded by the dashed box in Figure 1 are interconnected. In fact it has 
been shown that the links between diffusion and energetics are extremely im- 
portant in the seasonal variations of the thermosphere where the diffusion 
process can account for an increase in the temperature amplitude by as much as 
a factor of three (Mayr and Volland, 1972b). 

The efficiency for the excitation of variations in the thermosphere depends 
on a number of properties such as those emphasized in Figure 2. These are 
related to 

(a) heat conductivity 

(b) ion collision frequency and viscosity 

(c) eddy and molecular diffusion coefficients and 

(d) horizontal, vertical and time scales of Q. 

The significance of these properties becomes apparent 
simple algebraic equations which describe the continuity of 
and energy. 

Attributing characteristic lengths to the energy input’s 
and horizontal (^) and vertical (h) scales, and replacing all 
by 

3 1 13 1 3 1 

3 t a r ' r 30 * l' 37 “ h 


when looking at 
momentum, mass 

time duration (t) 
differential operators 

( 1 ) 
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where 

t = time 
6 - 1 at i tude 


z = height, 

the continuity equations take the simplified forms 


ApocApacAT 


for the vertical momentum equation, 


Ua 



for the horizontal momentum equation 


A p 



for the continuity equation and 


( 2 ) 


(3) 


( 4 ) 


0 



for the energy equation. 


( 5 ) 


In the equations (2) through (5) all constants and quantities have been omit- 
ted except those essential to identify a particular process. So c characterizes 
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a constant proportional to the mass m of the constituent, v ion collisions, e 
viscosity, H scale height, K heat conductivity, k Boltzmann’s constant 
and i indicates a phase difference of 90° due to the time derivative. A p , AT and 
Ap are the relative variations of mass density, temperature and pressure re- 
spectively, p and p are the time average pressure and mass density and U and 
W are the meridional and vertical velocities respectively. Q is the amplitude of 
the heat input. The effects from the Coriolis force and thus zonal winds are not 
considered in this discussion. In the denominator of Eq. (5) the first term rep- 
resents heat conduction, the second one represents advection and the third term 
represents the energy contributions that go into breathing and internal energy. 

For a minor constituent with a relative density amplitude Ap m and a vertical 
scale height H m the diffusion velocity W D = W m - W is derived from the continuity 
equation 


W « __! l H » h J i . T i 1 h „ 

”d ~ 


J_ _ J. 
A K 


1 _ J_ 

H h 

m m 


( 6 ) 


where h is the scale height of W , and its density is 


W 

A p «c AT-i 

m m j) 


( 7 ) 


where D is a factor proportional to the diffusion coefficient. Combinding Eq. (6) 
with Eq. (7) yields 
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Assuming that in first approximation the diffusion process does not affect 
the total internal energy or the total pressure variations, we can estimate the 
diffusion effect upon the temperature. Let AT be the temperature variation 
under the assumption of diffusive equilibrium and AT D the temperature varia- 
tion with diffusion; then the relation 


A T n a A T + 


(c p + c p ) D 

v rrr m 1 m ' 


can be derived neglecting for simplicity the temperature components in the pres- 


sure variations. AT is due to the diffusive equilibrium assumption of course a 
fictitious quantity while at d is considered to be a realistic representation of the 


actual variations in the gas temperature T g . 


Although heat conduction is the primary energy transport mechanism for the 
spatially uniform component in the thermosphere, it plays a relatively minor role 
in spatically non uniform components where the thermospheric circulation com- 
petes in the redistribution of energy. Thus the second term in the denominator of 
Eq. 5, describing the advective energy loss, is generally larger than the first one 
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which describes the heat conduction loss. Defining a characteristic time for the 
advective energy exchange 



_P 

P 


Eq. 5 can be further simplified to yield for the latitude dependent component 


for r »t and 

C 


for r «t . 

C 



( 10 ) 


( 11 ) 


r c is generally of the order of half a day; it depends however upon the elec- 
tron density as well as on the scale height, h,and the thermospheric density, p. 
Since the time duration of magnetic storms is typically of the order one or two 
days, relation (10) is approximately applicable. This shows then that the ampli- 
tude of the latitude dependent component tends to be proportional to the square 
of the horizontal scale, Z , and to the ion collision frequency V. The character- 
istics of the density variations, expressed in (10) and (11) have been discussed 
in detail by Volland and Mayr (1971). 


Substituting Eq. 10 into (3) and (6) the following relations can be deduced 


for the wind field 
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Ua i H 




(12) 



(13) 


which showsthat the vertical wind velocity tends to be independent of the hori- 
zontal scale i . Similarly, the diffusion effect for the minor constituents (see 
Eq. 8), which is proportional to W,is then also independent of the horizontal scale, 
and the variations of the minor constituents can be expressed in the form 


A P oc c AT 

'm ^ m 



(14) 


Considering Eqs. (2) and (10) the temperature variation, including diffusion 
(see Eq. 9), takesthe final form 

r , r\ 


A T oo J 

D 'v A 


■e 2 1 A .I 

H h 


V + 


p 




Pn 


hV P D (co+c PJ 


H H h + h 

^ m mj 

(tx) 


p 


(15) 


V J 

The first term on the righthand side of Eq. (15) represents the "tempera- 
ture" variation that would result under the assumption of diffusive equilibrium. 
The second term represents the diffusion effect which is independent of i and 
is relatively more significant for small scale disturbances where the amplitude 
of the "diffusive equilibrium temperature" becomes small. 
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Although small scale structures cannot be maintained in the total mass density 
(Eq. 10) they prevail in the vertical wind velocity and consequently tend to 
maintain the small scale structures in the individual constituents (second 
term on the righthand side of Eq. 14) and in the temperature (Eq. 15). 

The Equations (10) through (15) are strictly speaking only valid for periods 
in the disturbance which are much longer than the characteristic times for 
diffusion and advective cooling. For shorter periods the ^ dependences of course 
disappear and transport processes become less effective, a characteristic that 
is very significant considering that the duration of storms is variable and that 
each storm constitutes an entire spectrum of frequencies. 

III. THEORETICAL MODEL 

The theoretical model used in our calculations has been described in Mayr 
and Volland (1972b). It is essentially a three-dimensional, two-component model 
in which heat conduction, advective energy exchange, ion drag, viscosity and diffusion 
are considered in a self consistent form. Two versions of the two-component 
model are used. One in which the O and N 2 constituents are considered and 
which also provides a realistic description of the temperature, and a second 
one in which He is considered diffusing through a fictitious constituent with a 
mass corresponding to the mean molecular mass of all the other species. 

This model provides also an estimate of the "diffusive equilibrium temperature." 
To simplify the mathematical analysis perturbation theory is applied and the 
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latitude and time dependences are described in terms of spherical and Fourier 
harmonics respectively. 

In our calculations we adopted the same input parameters that have been 
used in Mayr and Volland (1972b) except for the height distribution of the heat 
input. Cole (1962, 1963, 1971) has shown that Joule heating is an important heat- 
ing mechanism during magnetic storms with the maximum of the heat input rate 
occuring at altitudes as high as 150 km. We have chosen therefore in our cal- 
culations a heat input distribution shown in Figure 3. 

IV. MAGNETIC STORM ANALYSIS 

To illustrate and verify the effects that have been qualitatively discussed 
in Section II we have computed the density, temperature and wind fields for 
heat inputs with various latitude structures described in the forms of the 
sperical harmonics P 0 , P 2 » P 4 and P 8 * The magnitude of the heat input at the 
poles is given in Figure 3. The time period chosen for these calculations was two 
days representative of the time scale for the major energy input of a magnetic storm. 

Figures 4 and 5 describe the relative amplitudes and phases for the three 
constituents N 2 , O and He and for the temperatures T g as well as the mass 
density respectively. In both figures the density scales are the same but the 
temperature scales are expanded by a factor of twenty. 

The distributions for N, , O and T have been derived from the N -O model, 

z g 2 

and since N 2 and O are the major constituents in the thermosphere the 


b 
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temperature variations from this model are considered to be realistic. The He 
distribution and the "diffusive equilibrium temperature" were derived from the 
He model in which the composition of the major constituents was assumed to be 
constant with time, n are the orders of the spherical harmonics P n and since with 
increasing n number the structure in the latitudinal variation increases it is ap- 
proximately related to the horizontal scale factor -t (see Section H) by 

1 

na I 

It is apparent from Figure 4 that the amplitude in N 2 decreases with in- 
creasing n in accord with the temperature decrease in Figure 5. The most 
drastic decrease occurs between P Q and P 2 thus demonstrating the importance 
of theadvective energy exchange associated with the thermospheric circulation 
set up in P 2 but not in P Q . Furthermore it can also be seen that the rate of 
this decrease slows significantly down at higher wave numbers n. This latter 
characteristic was expressed in Eq. 15! The "diffusive equilibrium contribu- 
tion" in the temperature decreases with increasing n (or decreasing horizontal 
scale -£) but the diffusion effect upon temperature and major constituent remains 
essentially independent of n (for n > 0) and thus becomes dominant at greater 
wave numbers. 

The importance of the diffusion effect for the minor constituents O and He 
is evident in Figure 4. In the globally uniform component P Q the diffusion 
process can only be excited by the one dimensional thermal expansion and 
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contraction of the thermosphere. Combined with the relatively large tempera- 
ture effect, diffusion is thus rather unimportant. This is apparent in both con- 
stituents O and He (solid lines) which vary more or less in phase with the 
temperature above their isopycnic levels at 175 and 300 km respectively with 
their amplitudes being essentially proportional to their mass at exospheric 
heights. 

For the nonuniform components P 2 through P 8 the thermospheric circula- 
tion is significantly involved in the diffusive redistribution of the minor con- 
stituents which, combined with the relatively small temperature effects, becomes 
very important. Thus both constituents are by almost 24 hours (180°) out of 
phase with the temperature and the N 2 concentration, implying that the minor 
species decrease at times when N„ and T increase, in agreement with our 

1 g 

earlier results (Mayr and Volland, 1972a). It is also apparent that this diffusion 
effect is essentially independent of n and thus contributes to the above discussed 
diffusion effect in the temperature. 

There are, however, some significant differences between the He and O 
distributions which bear upon the differences in the mass and abundance of both 
species. Accordingly, the diffusion effect is much less pronounced for O than 
for He. In addition the temperature effect is more important for O than for He 
and thus compensates the diffusion effect in O at higher altitudes. The conse- 
quence of the latter is that the amplitude in O initially decreases above 200 km 
with its phase gradually shifting toward the temperature phase at greater heights. 
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After going through an isopycnic level at exospheric heights the amplitude in 
O again increases with height and this constituent begins to vary, although greatly 
damped, more or less in phase with the temperature and the N 2 concentration. 
The height of this isopicnic level depends on the temperature amplitude. Thus 
it is lowest for P 2 (300 km) (and still lower of course for P Q (175 km)) and it in- 
creases to altitudes above 450 km for P g . For He with the weaker temperature 
and stronger diffusion-effect such an inversion does not take place (except for 
P 0 ) and its 180° phase difference with N 2 and T g remains intact throughout the 
thermosphere. 

It is interesting to note the rather complex behavior in the phase relations 
between the various physical parameters. As n increases the phases of the mass 
density and "diffusive equilibrium temperature" are advanced to earlier times at 
exospheric heights. This phase shift is particularly large between P Q and P 2 
when for n > 0 advection becomes very important thereby reducing 
substantially the characteristic time for the energy exchange. This latter 

mechanism is also the reason for the initial phase shift between the P and P 

0 2 

components in the temperature and the N 2 concentration. With increasing wave 
numbers n,the phases in atomic oxygen and helium are also shifted toward 

earlier times but for n > 0 the N and T phases are advancing toward later 

2 g 

times. With increasing n the diffusion process leads to an increasing separation 
of the phases of the individual constituents away from the phase of the mass density 
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and "diffusive equilibrium temperature." Thereby, the heavier constituent and 
the temperature go into one and the lighter constituents go into the other direction , 
indicating a sort of diffusive separation in time space. 

Finally, it is shown in Figure 5 that between P Q and P g the amplitudes of 
the mass density and diffusive equilibrium temperature decrease by about a 
factor of 100 at exospheric heights in agreement with the n dependence in Eq. 10 
and with our earlier results from an analytical model (Volland and Mayr, 1971). 

It is therefore understandable that the satellite drag data, which measure the 
mass density, can only observe the global components in the magnetic storm 
response of the thermosphere quite contrary to the composition and temperature 
measurements which should be capable of "recognizing" the local character of 
the disturbance. 

Figure 6 shows the amplitudes and phases for the vertical and meridional 
transport velocities of the N 2 constituent. For n = 0, that is for the global com- 
ponent, the horizontal velocities are of course zero. However with increasing 
n the meridional velocity decreases and the vertical velocities remain 
essentially constant, both features being consistent with the relations (12) 
and (13). The latter characteristic is responsible for the fact that the 
diffusion effect is nearly independent of n as has been pointed out before. 

For n = 0 the vertical velocity is smaller than that for n > 0 in the lower ther- 
mosphere where diffusion is important, however, it is greater at higher altitudes 
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thus accounting for the larger temperature amplitude in P Q and for the lack of 
horizontal mass transport. It can be seen that the meridional transport veloci- 
ties change their direction at about 150 km blowing away from the heat source 
above this level and blowing toward the heat source below it with the center of the 
circulation apparently occuring near the height for the energy input maximum. 

The vertical transport velocities for atomic oxygen and helium are up to 
30% higher than those of the major constituents in the lower thermosphere. 

Since all the vertical transport velocities are almost in phase there, this im- 
plies that a downward drag force is applied with its maximum value near the 
time of the maximum energy input, and this is consistent with the depletion 
effect of the O and He during disturbed conditions. 

V. MAGNETIC STORM SYNTHESIS 

To visualize the magnetic strom response in the thermosphere we con- 
structed a heat source in the form 

3 8 

Q = Z] cos (na;t) L Qnm p * {9) 

n = 1 m=0 

for cl > = 2 7 t / t withr being the time period of two days, and for m = 0, 2, 4, 6, 8 
thus describing a heat source that is symmetrical with respect to the equator. 
Figure 7a shows the distribution of the heat input rate at the time of the input 
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maximum and Figure 7b shows the time variation of the heat source in the 
auroral zones. The shape for the height variation of the heat input is the same 
that was shown in Figure 3. The values in Figure 7 were taken from 140 km. 

Since we restrict ourselves to a rather limited number of harmonics the input 
distribution is not too realistic in particular with regard to the waves at low 
latitudes and before and after the storm, and they should be ignored in the 
interpretation of our results. 

In the following we shall present our results in the form of contour plots 
in time-latitude frames. To facilitate the comparison with the heat input distri- 
bution we show in Figure 8 again the heat source in a contour plot. The time 
scale is chosen such that the heat input maximum occurs at 00 UT. To em- 
phasize that some of the structures before and after the storm and at low 
latitudes are unrealistic we show them in thin lines. 

Figures 9a through 9e show at 200 and 450 km the distributions of the N 2 , 

O and He constituents as well as those of T and the mass density p respectively. 

g 

The contour densities were scaled with a factor s, normalized to 1 for N 2 at 450 
km, such that in each figure the number of contours is approximately constant. This 
serves the purposes to provide the optimum information on the temporal and 
latitudinal structure and to provide a means of assessing the magnitudes of the 
variations by comparing the scale factors s. For the actual magnitudes of the 
variations the numerical values associated with each contour have to be con- 
sidered, and in evaluating these numbers it has to be kept in mind that they 
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represent relative variations to be superimposed upon the time average and 
globally uniform component of the thermosphere. 

From Figures 9a through 9c the disparity in the magnetic storm response 
of the various atmospheric constituents becomes apparent. The N 2 concentra- 
tion increases within the auroral zone where most of the energy is deposited, 
and this behavior is there in accord with the temperature enhancement shown 
in Figure 9d. The time and latitude scales in the variations of N 2 , however, 
are substantially larger than those in the heat source (Fig. 8). This is 
particularly apparent from the lack of structure at low latitudes and after 
the storm, and from the elongation of the N 2 bulge toward the poles and 
toward later times. 

In contrast, O in Figure 9b is shown to decrease in the auroral zone at 
200 km with some trend to increase at low latitudes. While the latitudinal 
variations reflect the structure in the heat source, the time scale in the 
response of O is due to the long diffusion time even larger than that in N 2 . 
Completely different are the characteristics in the distribution of O at 450 km. 
Here during the storm atomic oxygen peaks at the poles and decreases gradually 
toward the equator resembling not at all the rather localized heat input into the 
auroral zone. The time scale at 450 km is at a given latitude, however, much 
narrower than at 200 km with the density maximum progressing to later times 
as one moves away from the auroral zone. A comparison between the storm 
time variations in N 2 and O at 450 km revealsthat N 2 increases by about a 
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factor of 4 between equator and auroral zones while O increases only by about 
20% in this latitude range in substantial agreement with the observations of 
Taeusch et al., 1971 and with our earlier results (Mayr and Volland, 1972a). 

With regard to the phase relations it can be seen that in the auroral zone at 
450 km O peaks at about 2 00 UT followed by T g at 3 00 UT and N 2 at 3 30 UT, a 
behavior basically very similar to the response sequence in the diurnal varia- 
tions of the composition (Mayr and Volland 1972c). 

Helium which is even more affected by the diffusion process than O is also shown 
to decrease in the auroral zone. The amplitudes smaller than or equal to 
-1 are of course unrealistic thus reflecting upon the limitations in the perturba- 
tion approach. However, apart from that our He model is also for other reasons 
rather crude since thermal diffusion has been neglected and since by consider- 
ing only two components the temperature variations in the He model are un- * 
realistically small. A proper description of He would require a model in 
which the diffusion and energy transport of the three constituents He, O and N 2 
are considered in a self consistent form. Nevertheless, our results indicate 
that magnetic storms can induce a substantial reduction in He at times when the 
N 2 concentration is enhanced. Since the temperature is not very effective for 
He, it cannot cancel the diffusion effect (as for O) and hence the He decrease in 
the auroral zone and the enhancement at low latitudes remain intact throughout 


the thermosphere. 



22 


Figure 9e finally shows the response in the mass density p which is observed 
in the satellite drag data. Both at 200 and at 450 km the latitudinal variations 
in p during the storm are relatively insignificant and do not at all reflect the 
local heat input into the auroral zone. In particular.it is apparent that the 
global response in p is relatively significant (60%) when compared with the 40% 
increase between equator and pole during the storm at 450 km. This predomi- 
nance in the global component - which is even more pronounced at 200 km-con- 
firms thus our earlier conclusions (Volland and Mayr, 1971) and appears to be 
in substantial agreement with the satellite drag data. The maximum in the mass 
density at 450 km occurs near 2 00 UT at high latitudes and shifts toward 8 00 
UT at the equator, again in basic agreement both with Volland and Mayr (1971) 
and with the average time lag of 6°° hrs deduced from satellite drag data 
(Jacchia et al. 1967). 

Figure 10 illustrates the difference between the actual computed gas tem- 
perature and a sort of "effective temperature" which we deduced under the as- 
sumption of diffusive equilibrium from the variations in the mass density at 450 
km. The latter "temperature is equivalent to the one that has been inferred 
from the satellite drag data. The delay of six hours chosen for this comparison 
corresponds approximately to the response time of the global component in the 
mass density. For Figure 10 the relative temperature variations were converted 
into absolute values with a pre-storm temperature of 1000°K. It is then ap- 
parent that the gas temperature increases only slightly at the equator but 
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increases by about 250°K from the equator toward the auroral zone where a 
pronounced peak develops, hi contrast the "effective temperature" increases by 
120°K at the equator reflecting upon a relatively strong global component and 
increases by only 40°K from the equator toward the pole where an elongated 
"temperature" bulge develops. While the magnetic storm response in the mass 
density can thus be reconciled with the variations in the composition as observed 
on OGO-6 the temperature derivation from the satellite drag data can definitely 
not be justified. 

Figure 11 shows contour plots for the 0/N 2 ratio at 200 km deduced 
from our model. This ratio is an important parameter for the magnetic storm 
response of the ionosphere since the electron density at and below the F 2 peak 
varies proportional to the ionization rate, which depends on O, and inversely 
proportional to the ion-molecule charge transfer rates, which depend on N 2 
and 0 2 . Our results show that the electron density in the F 2 region could 
decrease up to about a factor of three at high and mid latitudes and slightly 
increase at low latitudes, a result which verifies the hypothesis of Seaton 
(1956), Duncan (1969), Chandra and Herman (1969), and Obayashi and Matuura 
(1972), who postulated this kind of composition changes in interpreting the 
magnetic storm response of the ionosphere. Obayashi and Matuura (1972) also 
suggested that the composition effect could be induced by the thermospheric 
circulation. However, in addition to these composition changes which influence 
the ionpsphere through the ion chemistry, at least three other features in the 
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magnetic storm response of the neutral atmosphere could be of importance for 
the ionosphere. These are (a) the nearly global increase in the thermospheric 
density which will decrease the ionospheric diffusion coefficient and thus enhance 
the electron density and height of the F 2 -maximum (b) the effects of the changing 
composition and temperature upon the energetics of the ionosphere which are 
probably rather complex when considered in conjunction with the processes of 
non local heating and (c) the effects of thermospheric winds (Kohl and King (1967), 
Obayashi (1972)) which obtain during storms an additional equatorward component, 
shown in Figure 12, that would tend to decrease the electron density below the F 2 
maximum and to enhance it at and above the F 2 peak. 

Finally, Figures 12 and 13 show the meridional and vertical components of 
the velocity field for atomic oxygen. Our results show that below 250 km, where 
oxygen is the minor constituent, the differences for the vertical transport 
velocities of O and N 2 are less than 30% while the horizontal transport velocities 
for both constituents are essentially identical. For this reason we consider the 
velocities in Figures 12 and 13 to be representative of the wind field. From 
these figures it is then apparent that during disturbed conditions a large scale 
circulation is set up with air rising in auroral latitudes and descending at low 
latitudes. The meridional velocities attain thereby values as high as 150 m/sec 
at mid latitudes, sufficient to affect significantly the F 2 region ionization. The 
meridional velocity is shown in particular to increase gradually toward mid 
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latitudes thus reflection the predominance of the large scale components (see 
the ^ dependence in Eq. 12). hi contrast, the vertical velocity component tends 
to be independent of ^ (gee Eq. 13) and hence reveal the details in the heat input 
distribution to a much greater extent, the latter being also the reason for the 
abrupt changes in the neutral composition as discussed before. 

VI. SUMMARY AND CONCLUSION 

Considering simple dimensional relations, some of the most important 
characteristics in the magnetic storm response of the thermosphere were 
qualitatively discussed thereby showing that the diffusive redistribution of 
atomic oxygen can have effects upon the temperature structure. A 
quantitative three-dimensional model of the thermosphere dynamics was then 
presented. In this model energy and diffusive mass transport associated with 
the magnetic storm-induced thermospheric circulation are considered in a 

self-consistent form. Our results show that these processes indeed play a 
major role in the thermosphere and thus account for a number of composition, 

temperature and mass density phenomena which are observed in the satellite 
drag data and in the neutral composition measurements on OGO-6. In particular 
the model describes 

(a) the relatively steep temperature enhancement toward high latitudes 


during disturbed conditions , 
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(b) the depletion effects in atomic oxygen and helium coincident with the 
temperature increase , 

(c) the negligibly small latitudinal variation in atomic oxygen in the pres- 
ence of a significant increase in N 2 and a similarly drastic decrease in He 
toward the auroral zones at 450 km, 

(d) the gradual increase in the mass density p between equator and pole 
superimposed upon a relatively strong global enhancement in p during disturbed 
conditions, 

(e) the gradual phase shift in the response time of p from 2 to 8 hrs be- 
tween pole and equator and 

(f) the phase sequence between the individual atmospheric parameters with 
the maximum in O occuring closest to the time of the energy deposition followed 
by the maximum in T and that of N in the auroral zone at 450 km. 

g 2 

In producing a temperature increase of about 400°K which corresponds to 
the storm time condition specifically described here, an energy input rate com- 
parable to that due to the EUV was required. We adopted thereby a height dis- 
tribution for the heat input which we considered to be representative of Joule 
heating (Cole, 1971). 

Although it appears that our model can describe some of the major features 
in the magnetic storm response of the thermosphere ,it must still be qualified 
as a first order approach toward a self-consistent treatment of the diffusion and 
energy transport mechanisms in the thermosphere dynamics. The application 
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of perturbation theory and the restriction to the first few harmonics which were 
considered independently characterize our approximations. Furthermore we 
restricted ourselves to a two component model which is not entirely adequate 
for He. Apart from these limitations .which were imposed for mathematical con- 
venience, there are a number of processes which have not yet been considered. 
Ion drifts (associated with the magnetic storm induced electric fields) which 
couple momentum into the neutral atmosphere (Cole, 1971), the feedbacks from 
the magnetic storm response in the ionosphere and the effects of particle pre- 
cipitation which deposit their energy presumably within the lower thermosphere 
and mesosphere (Rees 1972) are just a few of the mechanisms to be further 
investigated. 

VII. REFERENCES 

Chandra, S. and J. R. Herman, F region ionization and heating during magnetic 
storms, Planet. Space Sci., JF7, 841, 1969. 

Cole, K. D., Joule heating of the upper atmosphere, Aust. J. Phys., 15, 223, 1962. 
Cole, K. D., Joule heating of the ionosphere over Halley Bay, Nature, Lond. 199 , 
444, 1963. 

Cole, K. D., Electrodynamic heating and movement of the thermosphere, Planet. 
Space Sci., JL9, 59, 1971. 

Duncan, R. A., F region seasonal and magnetic storm behavior, J. Atm. Terr. 
Phys., 31, 59, 1969. 


V 



28 


Jacchia, L. G., Two atmospheric effects in the orbital acceleration of artificial 
satellites, Nature, 183 , 526, 1959. 

Jacchia, L. G., J. Slowey and F. Veriani, Geomagnetic perturbations and 
upper atmospheric heating, J. Geophys. Res., 72, 1423, 1967. 

Kohl, H. and J. W. King, Atmospheric winds between 100 and 400 km and their 
effects on the ionosphere, J. Atm. Terr. Phys., 29, 1045, 1967. 

Mayr, H. G. and H. Volland, Temporal variations in the thermospheric composi- 
tion, E05, 51, 789, 1970. 

Mayr, H. G. and H. Volland, Magnetic storm effects in the neutral composition, 
Planet, Space Sci., 20, 379, 1972a. 

Mayr, H. G. and H. Volland, Theoretical model for the latitude dependence of 
the thermospheric annual and semiannual variations, NASA- Document 
X-621-72-235, J. Geophys. Res., in press 1972b. 

Mayr, H. G. and H. Volland, Two component model of the diurnal variations in the 
thermospheric composition, NASA Document X-621-72-233, J. Atm. Terr. 
Phys., in press 1972c. 

Obayashi, T. and N. Matuura, Theoretical model of F-region storms, Proc. 

Solar Terr. Phys., 4, 199, 1972. 

Obayashi, T., World wide electron density changes and associated thermo- 
spheric winds during an ionospheric storm, Planet. Space Sci., 20 , 511, 


1972. 



29 


Rees, M. M., private communication, 1972. 

Roemer, M., Geomagnetic activity effect derived from Explorer 9 data, Phil. 
Trans. Roy Soc. London, A, 262, 184, 1967. 

Seaton, M. J., A possible explanation of the drop in F-region critical densities 
accompanying major ionospheric storms, J. Atm. Terr. Phys., _8, 122, 

1956. 

Taeusch, D. R., G. R. Carignan and C. A. Reber, Neutral composition variation 
above 400 kilometers during a magnetic storm, J. Geophys., Res., 27 , 8318, 
1971. 

Volland, H and H. G. Mayr, Response of the Thermospheric density to auroral 
heating during geomagnetic disturbances, J. Geophys. Res., 76, 3764, 1971. 



30 


VHI. FIGURE CAPTIONS 

Figure 1. Schematic diagram for the physical processes that influence the 
temperature, composition and wind field of the thermosphere. 

Figure 2. Schematic diagram for the physical properties of the thermosphere 
which affect the temperature, composition and wind field. 

Figure 3. Height distribution for the heating rate considered to be representa- 
tive of Joule heating (Cole, 1971). 

Figure 4. Computed relative amplitudes and phases for N 2 , O and He at the 
poles associated with the spherical harmonics P 0 , P 2 , P 4 and P 8 . The 
heating rate distributions used in the calculations had the same latitude 
dependences with polar amplitudes equal to the input distribution in 
Figure 3. 

Figure 5. Same as in Figure 4 but for the mass density, T g and the "diffusive 
equilibrium temperature." The latter was derived under the assumption 
that the mean molecular mass is time independent which corresponds only 
approximately to the assumption of diffusive equilibrium. 

Figure 6. The same as in Figure 4 but for U and W. While the form in the 
latitudinal variations of W corresponds again to P , P 2 , P^ , and P g , those 
for U are different. They are in particular sin q P , sin q (P 3 + 3/7 P x ) 
and sin Q (P + 11/15 P g + 7/15 P 3 + 3/15 P } ) corresponding to the heat 
input distributions with P 2 , P 4 and P 8 respectively. 
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Figure 7. Global energy input distribution at 140 km. Figure 7a shows the lati- 
tude dependence at 00 UT when the maximum occurs. Figure 7b shows the 
temporal variations of the heat input in the auroral zone. The wave struc- 
tures before and after the storm and at low latitudes are of no significance. 
They are merely a manifestation of the limited numbers of harmonics 
used in the magnetic storm synthesis. 

Figure 8. Contour plots of the energy distribution at 140 km for the magnetic 
storm disturbance. 

Figures 9a through 9e. Contour plots of the relative variations in the magnetic 
storm response of N ,0, He, T and p calculated for 200 and 450 km with 

2 g 

the energy input distribution defined in Figures 3, 7 and 8. 

Figure 10. Computed gas temperature (solid line) in comparison with an "effec- 
tive temperature" (dashed line) deduced from the mass density under the 
assumption of diffusive equilibrium. 

Figure 11. Contour plots of the relative variations in the magnetic storm re- 
sponse of the 0/N 2 ratio of 200 km a parameter which is significant for the 
ionosphere. 

Figure 12. Contour plots of the meridional component in the transport velocity 
of O. Since the difference between the transport velocities of O and N 2 
are negligible below 250 km this velocity is considered to represent the 
wind velocity. 
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Figure 13. Contour plots of the vertical component in the transport velocity of 
O. The difference with the corresponding component in N 2 can be 30% and 
thus the values in this figure represent only approximately the wind 
component. 
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